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We use the generic replica symmetric cubic field-theory to study the transition of short range 
Ising spin glasses in a magnetic field around the upper critical dimension, d — 6. A novel fixed-point 
is found, in addition to the well-known zero magnetic field fixed-point, from the application of the 
renormalization group. In the spin glass limit, n — » 0, this fixed-point governs the critical behaviour 
of a class of systems characterised by a single cubic interaction parameter. For this universality 
class, the spin glass susceptibility diverges at criticality, whereas the longitudinal mode remains 
massive. The third mode, the so-called anomalous one, however, behaves unusually, having a jump 
at criticality. The physical consequences of this unusual behaviour are discussed, and a comparison 
with the conventional de Almeida- Thouless scenario presented. 
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The mean-field theory of the Ising spin glass JlJ pro- 
vided an astounding complexity of equilibrium prop- 
erties, showing up how disorder and frustration may 
lead to an unusual thermodynamics. More than two 
decades have elapsed since the ultrametric solution of 
the Sherrington-Kirkpatrick (2| model by Parisi was pub- 
lished in a series of papers (see jjj for the references), 
nevertheless consensus has not been reached about the 
validity of the mean-field picture for finite-dimensional, 
short range systems. The alternative scenario, the so- 
called droplet picture || ||, claims that the complex 
phase space structure is an artifact of mean-field the- 
ory, the glassy state consists of two phases related by the 
global inversion symmetry of the spins. 

The investigation of the spin glass transition in an ex- 
ternal magnetic field may resolve the debate: The glassy 
transition along the de Almeida-Thouless (AT) line || is 
a distinctive feature of mean-field theory, whereas spin 
glass ordering is destroyed by any nonzero magnetic field 
in the droplet model. Although a lot of numerical work 
has been performed, no convincing evidence has emerged 
until now in favour of either theory. An AT line was 
found in the four-dimensional case in Refs. j(| whereas 
numerical results in three dimensions were interpreted, 
although less convincingly, to support mean-field-like be- 
havior in || ||, |). On the other side, Ref. |l(J in- 
terprets the simulation data of j^j as quite consistent 
with droplet theory, and an analysis of the ground states 
in jllj showed that the spin glass phase of the three- 
dimensional model does not survive in any finite mag- 
netic field. More recently, however, an extensive study of 
the energy landscape Q suggests that a nonzero critical 
field may exist at zero temperature, separating the spin 



glass and paramagnetic phases. 

In this letter, replica field theory, as an alternative to 
numerical calculations, is used to attack the problem by 
extending the renormalization group study of Ref. p3[ . 
Our starting point is the Edwards- Anderson (EA) p4] 
model of N Ising spins on a d-dimensional hypercubic 
lattice, defined by the Hamiltonian 

H = -^JijSiSj - H^Si (1) 

(ij) » 

where the first (second) summation is over all nearest 
neighbour pairs (ij) (all lattice sites i), respectively. Jij 
are independent, Gaussian distributed random variables 
with mean zero and variance A 2 , and a homogeneous 
magnetic field H has also been included. The application 
of the replica trick followed by a Hubbard-Stratonovich 
transformation produces a replica symmetric field theory, 
with the Lagrangean £ m i cr taking over the role of TL. In 
this replica field theory the fields depend on two replica 
indices with the restriction <p a ^ — <j)^ a and (f) aa = 0; 
hence we have n(n — l)/2 field components. The replica 
number n must go to zero to reproduce quenched aver- 
ages; we will argue, however, that it is necessary to keep 
it finite until the very end of the calculations Q| . As it 
is common in the theory of phase transitions, the micro- 
scopic Lagrangean £ m i cr is replaced by an effective one 
C = CS 1 ^ + £( 3 \ obtained by iterating the renormaliza- 
tion group until irrelevant operators can be neglected. 
As a result of the replica trick, a generic replica symmet- 
ric field theory follows, which can be best represented 
in terms of operators invariant under the permutation of 
the n replicas: 
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The above field theoretical model is general enough to 
describe a large variety of possible transitions from the 
high temperature, replica symmetric phase, just below 
the upper critical dimension d = 6. The failure of Bray 
and Roberts Jl3| to detect a fixed point corresponding 
to the conventional AT transition gave definite support, 
from the analytic part, to the droplet theory. In this 
letter, we present a novel fixed point characterized by a 
rather unusual property, and propose it as the relevant 
one for the generic replica symmetric phase, i.e. that with 
a nonzero order parameter. We also put forward a possi- 
ble physical scenario whose validity is a prerequisite for 
that fixed point to control the spin glass transition in 
a magnetic field. Checking this, however, is out of the 
scope of this letter, and we leave it for future work. 

Even a leading order renormalization group calculation 
is blocked by the difficulties arising from the numerous 
and complicated replica summations, and from the fact 
that £( 2 ) is not in a diagonalised form. In Ref. |15j we 
worked out a transformation to a new set of bare pa- 
rameters (rn, ta and rj, for the masses, and g%, .. . ,gg 
for the cubic couplings) rendering the one-loop calcula- 
tion feasible. As an illustration, we computed in pUj ] 
the true masses T^, Ta and to one-loop order. To 
obtain the renormalized cubic interaction, we calculated, 
via long but relatively straightforward algebra, the trian- 
gle graph. In order to be completely parallel to Ref. jl3) , 
we chose the same renormalization scheme of integrating 
out degrees of freedom in the infinitesimal momentum 
shell between e~ dl A and A, where A is the ultraviolet 
cutoff. It is obvious that there is no sufficient space to 
present here the recursion relations in their total gener- 
ality, as worked out in Jig . The structure of these RG 
equations is as follows: 



dn/dl = TZi(r R ,r A ,r L ;gi, . 
dgi/dl = Gi(r R ,r A ,r L ;gi, . 



,g s ) i = R, AotL, 
g 8 ) i = l,..., 8. (4) 

The above set of equations must comprise two special 
cases, known from the literature for some time, providing 
us with a good check: 

1 . The zero magnetic field case was studied in Ref. |l7j] 
up to 0(e 3 ), the Lagrangean corresponding to it 



(fR = r A = tl = r, w\ = w and u>i — for i = 
2, . . . , 8) proves to be an invariant subspace of the 
set of Eqs. (U), expressing the higher symmetry this 
system possesses. The fixed point attracts a critical 
line inw-r space, which is totally massless (r^ = 
= Tl = 0), and, following the ideas in [fI5[ , 
we can identify this theory as the relevant field- 
theoretical model for the spin glass transition from 
the paramagnet to the generic replica symmetric 
phase with a nonzero order parameter. The fixed- 
point value w* 2 = e, with e = 6 — d, and 
the exponents 77 and v are in complete agreement, 
up to first order in e (see fig), with the results of 
Green for n = 0. As in mean-field theory, n is a 
rather innocent parameter around this fixed point. 

2. The equations of Bray and Roberts [|l3| are re- 
produced by assuming tentatively a critical surface 
with Tr = 0, whereas Ta and Tl finite. At some 
hypothetical fixed-point, the bare anomalous and 
longitudinal masses are infinite, resulting in a pair 
of recursion relations for the replicon-like couplings 
gi and g 2 : 



dgi/dl = g i {gi,g 2 ) 



= 1,2. 



(5) 



Using the relations gi = w\ and g 2 = "2w 2 [P-5| , we 
arrive at the RG recursions whose physically rele- 
vant fixed-point was searched for in vain by Bray 
and Roberts [^3| . 

A new theory, invariant under renormalization to O(e) 
order, emerges if the condition of degeneracy between 
the longitudinal (L) and anomalous (A) modes is re- 
moved. The physical picture behind this may be the 
following: The masses characterizing a critical mani- 
fold are observable through correlation functions. For 
the original lattice system of Eq. (^), three distinct spin 
glass correlation functions can be defined (see e.g. in 
||), their zero momentum limits will be denoted by G\, 
G 2 and G3. The transition in a field can be character- 
ized by the divergence of the spin glass susceptibility, 
Xsg = G\ — 2G 2 + G3 = T^ 1 , manifesting itself in the 
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criticality of the replicon (R) mode, while the longitudi- 
nal mode, G\ — 4G 2 + 3G3 = T^ 1 , remains massive. As 
for the anomalous one, we have the exact relation for n 
small but finite: 



l-nT L [G 2 - §G 3 ) 



(6) 



In mean-field theory, the combination G2 — IG3 is also 
analytical along the AT-line, leading to the following ra- 
tios for the leading singularities: Gi : G2 : G3 = 1 : 
1/2 : 1/3. Dimension-dependent subleading singularities 
of the form ~ occur in finite dimensions, t being 
the reduced temperature, and an analysis of the pertur- 
bation expansion of the propagators shows that for the 
next-to-leading term, \x is equal to 4 — d/2, d > 6. As 
a result, a jump of the anomalous mass may develop be- 
low d — 8 |p2| , provided the longitudinal mode remains 
massive. From Eq. (||) it follows: 



if n — > first, 
if t ->■ first. 



(7) 



From the assumptions Tl finite and /1 > 0, and using 
3), a jump in the anomalous mass necessarily follows 



in the limit n — ► even below d — 6. Whether or not 
these assumptions are correct, can be tested by the RG 
equations (|]). To detect this behavior, we now search for 
a nontrivial fixed point with Tr = T A = and = 00. 
It is obvious from Eq. (0) that n is a crucial parameter 
now, and it must be finite when computing the RG flows, 
setting it to only at the very end of the calculation p3| . 
The longitudinal bare mass is kept at its infinite fixed- 
point value, longitudinal- like couplings (174, 177 and g%; see 
Eq. (48) of Ref. |J|) decouple from the rest of Eq. @. A 
massless renormalization scheme (massless with respect 
to the replicon (R) and anomalous (A) masses) can be 
deduced from the remaining part of (pj): 

dgi/dl = Gi(gi,g 2) g3,g5,g6) i = 1,2,3,5,6. (8) 

In our one-loop calculation, Qn is a cubic polynomial of 
its variables with coefficients which are rational functions 
of n. It turns out that Q3 and Qq are always zero when- 
ever <?3 = <?6 = Jl6|, indicating that the 3-dimensional 
manifold so defined is an invariant subspace of the RG 
equations (Js[) . Instead of presenting long and cumber- 
some formulae for generic n, we display the recursion 
relations for gi, g 2 and 175 only in the n — * limit: 
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(r]A above is proportional to n, i.e. it is zero here. It is 
shown only to display the generic structure of the RG 
equations.) 

We found a novel non-trivial fixed-point from Eqs. 
(§), (g) and (§): 



fJi 



e/2, 



g ; = -v^/4. (10) 



The existence of this fixed-point is due to the term and 
g 2 in Eqs. ( |9a| ) and (|9q), respectively; omitting them, we 
just get back the Bray-Roberts equations (||). It is re- 
markable that, in a sense, we have found in Eqs. (9a), 
(ph|), ( pig 1 ) and (9d) a generalization of the RG theory 
put forward in Ref. p3[ . We must notice, however, that 
beside the replicon mode, the anomalous one is also crit- 
ical on the manifold attracted by the fixed-point (p"o|). 
The most striking feature of Eqs. (|9a|), ([)t]) and (|9c|) 



is that they coincide for g\ = w, g 2 = 2w and 95 = —w/2 
providing the single parameter RG equation for w: 



dw e 
-dl=2™ 
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Translating this to the language of the w couplings in Eq. 
(^), w\ = w 2 = wq = w, and all the other w's are zero. 
The system with the single cubic operator 
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rescales under RG, evolving into the fixed point w* = 
■s/e/2, where the corresponding eigenvalue is X a , = — e. 
This result can be compared with the zero-field case, 
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where the similar rescaling property of the system ex- 
pressed a higher symmetry than the permutation invari- 
ance of the n replicas. There is one important distinction 
we must notice, however: the rescaling behavior of the 
zero-field Lagrangean under iteration is independent of 
n, whereas it develops only in the spin glass limit n — > 
for the system with the cubic coupling in (p"l|). 

Including the masses into the RG scheme, it can be eas- 
ily checked that the condition tr = va = f is preserved 
under iteration. Using results from Ref. |t5), namely Eqs. 
(22)-(24) and (28)-(30), mi = f/2 and m 2 = follow 
then, while 7713 is infinite, inducing the freezing-out of the 
longitudinal component of <f a ^ . The quadratic operator 
in the brackets of Eq. || reduces to the simple replicon-like 
invariant 
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although <fi a P has now an anomalous component too. 

We can deduce critical indices belonging to this new 
fixed point, and we display them here for completeness: 
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To connect these to usual exponents like that of the spin 
glass susceptibility, \sg ~ i~ 7 , is not trivial now due to 
the coexistence of two critical masses, and needs further 
study. 

We propose the simple model of Eqs. ([l2|) and ( fTi"| ) as 
a candidate for studying the replica symmetry breaking 
(RSB) transition from the replica symmetric phase with 
a nonzero order parameter; the spin glass transition in 
an external field belongs to this class. Nevertheless, it is 
difficult to find evidence for this. The check by testing 
the crossover, from a £ m icr in the vicinity of the zero- 
field critical point to the new fixed point, is blocked by 
the large distance in this huge parameter space, and by 
the evolution of the longitudinal mass from the near-zero 
value to infinity. It is also obvious that the irrelevant 
operators present in £ m icr influence this crossover, ren- 
dering this check very difficult. To by-pass this problem, 
it is tempting to imagine an alternative scenario, viz. the 
existence of the replica symmetric phase with nonzero 
order parameter even in zero field. This two-step pro- 
cess from the paramagnet to the RSB phase is present in 
mean-field, but only for finite, albeit infinitesimal, n [ p"5| . 
If this scenario occured in low enough dimensions even 
for n = 0, the crossover from one type of transition to 
the other would disappear. In this case, the intermediate 
replica symmetric phase with Q =/= would have the re- 
semblance to a droplet-like phase. It is clearly necessary 
to perform further investigations, mainly a higher order 
calculation and numerical investigations. 

In conclusion, we must stress that the theory we have 
put forward for the RSB transition is qualitatively dif- 



ferent from the AT transition of mean-field theory. We 
argued that the change occurs at d = 8, below which 
the relevant Gaussian theory is that with zero replicon 
and anomalous masses, while infinite longitudinal one. 
It is this Gaussian fixed-point which gives birth to the 
nontrivial one we found, governing the RSB transition 
below d — 6. We can speculate that these qualitative 
differences may affect the glassy phase too, resulting in 
a more general RSB scheme than the ultrametric one of 
mean-field theory. 
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grateful to Imre Kondor for a critical reading of the 
manuscript. This work has been supported by the Hun- 
garian Science Fund (OTKA), Grant No. T032424. 
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